Magnetoelectric effects in superconducting nanowires with Rashba spin-orbit coupling 
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Recent experiments in semiconductor nanowires with a spin-orbit coupling and proximity-induced 
superconductivity exhibit signatures of Majorana bound states predicted to exist in the topological 
phase. In this work we predict that these nanowire systems exhibit unconventional magnetoelec- 
tric effects showing a sharp crossover behavior at the topological phase transition. We find that 
magnetic fields with a component parallel to the spin-orbit field can give rise to currents in equi- 
librium. Surprisingly, also fields perpendicular to the spin-orbit field may induce currents and can 
be employed in adiabatic charge pumping. The perpendicular field magnetoelectric effect may be 
regarded as a manifestation of the anomalous Hall effect in one dimension. We discuss how the 
predicted phenomena could be observed in experiments and employed in probing the topological 
phase transition. 
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Introduction- Recent experiments in nanowires with a 
strong spin-orbit coupling and proximity superconductiv- 
ity show intriguing signatures of Majorana bound states 
[lH3j, a particle-like many-body excitations behaving as 
their own antiparticles. The surge of experimental and 
theoretical activity around the topic was largely initiated 
by the prediction that these systems exhibit topological 
superconductivity 0, [3 ■ Shortly before it was proposed 
that two-dimensional magnetic Rashba systems in the 
proximity of a superconductor provided a candidate to 
realize topological superconductivity [5, 6]. However, tai- 
loring suitable systems and controlling relevant physical 
parameters have proven more feasible in nanowires [7|- 
9]. Nanowire systems have also opened up exciting pos- 
sibilities to experimentally study effects resulting from 
the three necessary ingredients of the topological phase, 
a spin-orbit coupling, superconductivity and magnetic 
fields. 

In this work we predict and study properties of mag- 
netoelectric effects in Id wires resulting from the inter- 
play of the Rashba coupling, superconductivity and mag- 
netic fields. These magnetoelectric effects lead to non- 
dissipative electric currents that are not related to the 
phase gradient of the superconducting order parameter. 
The magnetoelectric effects exhibit a clear crossover be- 
havior when the nanowire undergoes a topological phase 
transition, thus providing a method to probe the tran- 
sition. We study two qualitatively different effects dis- 
tinguished by whether the applied magnetic field has a 
finite component parallel to the spin-orbit field or not. 
Experiments have demonstrated that the direction of the 
spin-orbit field can be resolved quite accurately 0| ■ 

First we predict that magnetic fields having a finite 
component parallel to the spin-orbit field lead to a finite 
DC current through the wire. This effect is sensitive to 
the closing of the energy gap associated to the topologi- 
cal phase transition. Then we show that electric currents 
can also be driven by magnetic fields perpendicular to the 




Figure 1: a): Superconducting nanowire with a Rashba 
spin-orbit coupling in external magnetic field B where 
Bso denotes the spin-orbit field for particles moving in 
the positive y direction, b): Experimental setup for 
probing the magnetoelectric effects. 



spin-orbit field. This phenomenon arises from the same 
band curvature effect that contributes to the anomalous 
Hall effect in 2d Rashba systems and enables adiabatic 
charge pumping when magnetic field is periodically mod- 
ulated. We discuss how the predicted magnetoelectric 
effects could be observed in experiments. 

Model system- We consider a nanowire with a Rashba 
spin-orbit coupling in the proximity of a superconductor 
in magnetic field, depicted in Fig. [T] The wire is de- 
scribed by the Bogoliubov-de Gennes Hamiltonian [|| , 



H{k) = (e k + ak<T x ) t z + B ■ a + At x 



(1) 



j-2.2 

where et = Aj— n and <Ji and Tj are Pauli matrices 

operating in the spin and the particle-hole space, respec- 
tively. Hamiltonian |lj is written in the Nambu basis 

* = (i'ktr^kii V'-fcn — i )1< -k^) T an d the coordinates have 
been chosen so that k is the momentum along the y-axis 
and the effective spin-orbit field is parallel to x direction. 
The first term corresponds to kinetic energy of electrons 
and holes including the Rashba coupling characterized by 
a coupling constant a. The second term is the Zeeman 
coupling B = (B x , By, B z ) due to the applied magnetic 
field and the last term arises from the proximity-induced 
superconducting pairing. Hamiltonian ([lj possesses a 
particle-hole symmetry and the four energy bands satisfy 
EL,$(k) = — Ei{— k) for i = 1,2 where the corresponding 
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eigenstates i and —i are related by a particle-hole trans- 
formation. In the case B x = 0, the spectrum admits an 
analytical solution 



Ei /2 QL)=4 + a*k*+B 



A 2 



T2jB 2 (e 2 +A 2 )+e 2 a 2 k 2 , (2) 



where B 2 = B 2 + B 2 . We adopt a convention accord- 
ing which B x (that is parallel to the effective spin-orbit 
field) is called parallel magnetic field and B y and B z are 
called perpendicular fields. For the reference, the rele- 
vant physical parameters in experiment employing InSb 
wires in Ref. [H were m — 0.015 m e , Er = 
AteV, A = 250 ^eV (at T=60 mK) and |B| 



^ = 200 
< 1.5 meV 

(for fields < IT). The transverse mode separation in the 
wire was several meV which is much larger than the other 
energy scales. 

Magnetoelectric effect for parallel fields B x ^ 0- In 
this section we consider magnetic fields that have a non- 
vanishing parallel component with the spin-orbit field 
(according to conventions of Fig. [1] this means that 
B x ^ 0). This leads to a non- vanishing current through 
the wire, an effect which is similar to the magnetoelectric 
effect predicted by Edelstein in 2d systems 1(1 XI |. 



Current through the wire can be evaluated as 



k.i 
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where J 
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dhHr z is the current 



operator, Pi is a projection operator to the eigenstate 
\rii) of Hamiltonian {T]) with energy Ei and /j is the 
Fermi distribution at energy Ei. The projection oper- 
ators are given by Pj = Ylj& g -g f° r tne f° ur bands 
and the trace should be calculated over the Nambu and 
spin indices. A further evaluation of Eq. © yields 
J = I / § E 4 + |Tr [P % a x ]) fi. The general formula 
for the current in terms of the energies is lengthy and 
presented in the supplement. Current is an odd function 
of the parallel field J(B X ) = — J(—B x ), a property which 
can be employed in distinguishing the magnetoelectric 
contribution. In the linear order in B x current is given 

by 
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A 2 +s 2 +a 2 k 2 + (-l) 



(4) 



where D = A y /B 2 {e 2 , + A 2 ) + e 2 z a 2 k 2 . The wire under- 
goes a topological phase transition when the perpendic- 
ular field satisfies B = y/ /i 2 + A 2 [7|, associated with 
the energy gap closing E±\(k) = at k = 0. The numer- 
ator in the integrand for i — 1 term also vanishes at that 
point and current remains analytic. 




Figure 2: a): Current from Eq. (j?)) a function of the 
perpendicular field B for T = 0.025 Er, (i = 0.0, and 
A = 0.8, 1.0, 1.2, 1.5 x E R (top to bottom), b): 
Differential current ^ calculated from the curves in a). 
The maxima coincide with the critical field above which 
the system is in the topological phase. 



Current Q is plotted MS c\ function of the perpendic- 
ular field B in Fig. [2] a). Current starts from zero at 
-6 = and increases rapidly until saturating and start- 
ing to decrease slowly. Due to the combined effect of the 
spin-orbit coupling and magnetic field B Xl the disper- 
sion becomes asymmetric Ei(k) ^ Ei{—k). This asym- 
metry and superconductivity are necessary ingredients to 
achieve finite current since in the normal state the current 
cancels [lj|. For large B, the eigenstates resemble mag- 
netic rather than helical metal. For a dominantly mag- 
netic dispersion B x induced k < — > — k spectral asymme- 
try becomes suppressed, thus explaining the saturation 
and decrease of current for large B. The saturation sets 
on when 8b J achieves maximum which, as illustrated 
in Fig. [2] b), coincides accurately with the topological 
phase transition. The results in Fig. [5] are calculated at 
temperature which translates to 60 mK for parameters of 
Ref. [l|. The peaks in Fig.[5]b) get broadened and shorter 
as the temperature increases, but the well-defined struc- 
ture survives to well above T = 0.1 Er. Thus the par- 
allel magnetoelectric effect provides a method to probe 
the phase transition indirectly. It is important to find 
alternative methods to characterize the phase transition 
since there exists theoretical evidence that some of the 
considered signatures of the topological phase, such as 
the zero bias peak in the tunneling current 12j and the 
fractional Josephson effect [13j| . could have explanations 
different from Majorana bound states. 

So far we have analyzed translation-invariant systems 
but experiments are performed with finite wires. For ex- 
perimental implications we need to consider a setup in 
Fig. [U b) where a finite wire is laterally coupled to two 
reservoirs L and R. Superconductivity in the wire is in- 
duced either by one or both of the reservoirs. Both of the 
reservoirs need to be superconducting in order to achieve 
finite current since equilibrium currents cannot be sus- 
tained in a normal metal reservoir. In practice the wire 
may support multiple transverse modes as in Ref. [lj. 
However, when the mode separation is large compared to 
the other energy scales, as in [l|, modes can be thought 
as independent (at least in the clean limit). In addition, 
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the modes below the topmost one do not contribute to 
current when the mode separation is much larger than 
other characteristic energy scales, so the single mode ap- 
proximation is expected to be an accurate starting point. 

It is also interesting to consider the case where one of 
the reservoirs is in the normal state or the wire is effec- 
tively pinched of by a local gate so that current is sup- 
pressed. According to thermodynamic arguments con- 
sidered in the 2d case [ 1 ljj , the magnetic field induced 
current is canceled by a superconducting phase gradient 
along the wire in equilibrium. Thus the effect could be 
perhaps also observed by f^-dependent build up of su- 
perconducting phase when current is blocked. 

Magnetoelectric effect for perpendicular fields: linear 
response - Now we consider magnetic fields that are per- 
pendicular to the spin-orbit field which in our conventions 
mean that B lies in the y — z plane (B x — 0). Below we 
derive the surprising result that perpendicular fields can 
also be employed to drive currets in the wire. This mag- 
netoelectric effect arises from essentially the same band 
curvature effect that gives rise to the anomalous Hall 
effect in 2d Rashba systems with perpendicular magne- 
tization pi. 

First we consider a magnetic field configuration with 
a static component B in z-direction and a weak time- 
dependent component SB y in y-direction. Treating SB y 
as a perturbation, current in the wire is given by the 
standard linear response theory as 

J(u)± = x(u)SB v (u). (5) 

The imaginary-time representation of the response func- 
tion x{ u ) i s given by 

X(ivm) = jo ^ Tr[JG(i(u m + Vn))<7yG{iv n )], (6) 

" k,n 

where (3 is inverse temperature and the Matsubara 
Green's functions G{iv n ) are defined by G{iv n ) = . X _ H , 
where H is given by Eq. JT]) with B = (0, 0, B). Expres- 
sion (j6|) can be evaluated using methods of Ref. [H} where 
similar expressions were evaluated in 2d geometry. In the 
following we are interested in low-frequency dynamics. In 
lu — > limit x{ u ) approaches to zero so it is convenient 
to consider the quantity lim LJ _ ! .o x(cj)/icj = — J T^£l xy , 
where 

v ' n,n f 

Here P n is a projection operator to the eigenstate with 
eigenvalue E n given by Eq. ([2]). The explicit form of 
n xy (k) is long and given in the supplement. Current in 
the wire can now be expressed as 

'«>-^(!/§«w»>)- m 




Figure 3: a): Magnitude of the perpendicular effect as 
a function of the field B for T = 0.025 E R , n = 0.0 and 
A = 0.8, 1.0, 1.2, 1.5 x E R (from top to bottom), b): 
Adiabatic pumping scheme. When the perpendicular 
field B traverses a loop in the y — z plane, charge Q is 
pushed through a cross section. 



Thus we have reached a remarkable conclusion that tem- 
poral variations of magnetic field parallel to the wire can 
drive current. Figure [3] a) illustrates the behavior of di- 
mensionless quantity — f ^£l xy setting the magnitude of 
the effect. At low temperatures it behaves essentially as 
the current from the parallel magnetoelectric effect plot- 
ted in Fig. ||2j a), therefore the differential current 8b J 
resulting from the perpendicular effect (jHJ exhibits peaks 
similar to those shown in Fig. @ b). From symmetry 
considerations it follows that result ([5]) holds even if the 
directions of the static and dynamic magnetic fields B 
and 8B y are rotated arbitrarily around the s-axis as far 
as they remain perpendicular to each other. We are not 
aware of previous predictions of constitutive relations be- 
tween current and magnetic field of type (|5]) in Id wires. 

The magnetoelectric effect (J8]) can be understood as 
a manifestation of the same band curvature effect that 
gives rise to the anomalous Hall effect in 2d Rashba sys- 
tems [lH • Recently the anomalous Hall conductivity was 
calculated in superconducting two-dimensional Rashba 
systems with magnetization perpendicular to plane (l5| . 
The Hall conductivity is given by o xy — \ j (2tt)^ ^xy(k) 
where momentum argument should be interpreted as the 
magnitude of the in-plane momentum k = ^Jk x + ky 
and the static part of the magnetic field B corresponds 
to magnetization perpendicular to the plane. Therefore 
&xy{k) can be regarded as the superconducting counter- 
part of the Berry curvature in non-superconducting sys- 
tem [Tl . \vi\ . The reason why Eq. (JSj) is connected to the 
2d Hall conductivity can be traced to the formal similar- 
ity of response function © and the off-diagonal current 
correlation function in 2d magnetic Rashba system [19j. 

The magnetoelectric effect (jSJ could be observed in a 
setup of Fig. Q] b) where both reservoirs are supercon- 
ducting. For example, by applying a static field and a 
small AC field perpendicular it, both lying in the y — z 
plane, the wire should exhibit a measurable AC current. 
Current as a function of the static field also exhibits a 
clear crossover behavior at the topological phase transi- 
tion and could be employed to probe it. 

Magnetoelectric effect for perpendicular fields: adia- 
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batic pumping- The relation between current and per- 
pendicular magnetic field JSJ suggests that the integrated 
charge flowing through a cross section when the field is 
modified depends only on the net change of the field but 
not on time that the procedure lasts. This is a hallmark 
of adiabatic transport and geometric pumping (l6l . 17 1. 



Below we discuss pumping processes where magnetic field 
executes a cycle in the y — z plane (B x = throughout 
this section), as illustrated in a particular case in Fig. [3] 
b). The calculation proceeds along the lines of adiabatic 
transport and polarization of insulators [17| . 

Suppose that some parameter A(t) in Hamiltonian (TJJ) 
varies slowly in time H(t) = H(\(t)). Let us further 
assume that the system is in the ground state and that 
the population of the energy band E n is f n (for negative 
energy bands /„ = 1 and for the positive energy bands 
In = 0) initially at t = 0. Current in the wire for t > 
is given by 



J (*) = / ^E^(*)i*(*)>/»> 



(9) 



where again J = ^dkHr z . Initially the states \ifj n (t)) 
satisfy \ip n (0)) = \n(0)) where |n(0)) is an eigenstates of 
(JTJ at t = 0. According to the adiabatic theorem, the 
temporal evolution of an eigenstate is given by \ip n (t)) = 



-if*E„{s)ds 



(W*)>-i&En'*. 



\n'(t))(n'(t)\d t n) 



the 



lowest order in the time derivatives. Here \n(t)} and E n 
correspond to the instantaneous eigenstates and eigenval- 
ues of the time-dependent Hamiltonian. Following the 
the standard steps in adiabatic transport theory, when 
the adiabatic parameter traces path from A(0) = A,; 
through X(Tf) = A/, the total charge Q 
pumped through the wire is 



fo f dtJ(t) 



Q 



Xi 



Tr [d k HT z P n ,d x HP n ] 
(E n -E n ,y 



(/„-/„-)• (io) 



Here we have introduced projection operators P n {X) to 
instantaneous eigenstates corresponding to E n {X). The 
integrand in Eq. (JTUJ) is analogous to the Berry curvature 
in non-superconducting systems [l7j. The extra Nambu 
matrix t z inside the trace, arising from the Nambu struc- 
ture of the current operator, makes the analogy to the 
Berry curvature incomplete and the pumped charge is 
not quantized. Result (fT0|) can be generalized to finite 
temperatures by replacing /„ by Fermi distribution at en- 
ergy E n (\ = 0) if energies E n do not change during the 
process or the cycle is faster than relaxation processes 
(while being slow compared to the relevant spectral gaps 
as the adiabatic approximation requires). 

Let us now consider the pumped charge when magnetic 
field is modified slowly in the y — z plane. The cycle 



executed by the magnetic field is conveniently parame- 
terized by magnitude B(X) and angle ^(A) so that field- 
dependent part of Eq. (p} becomes B ■ a — B cos(ip)a y + 
B sm((p) a z . Application of Eq. (fT0|) yields 



Q = 



d\Bd\tp 



dk 
2tt' 



^xy (&) ? 



(ii) 



were Sl xy {k) is given by Eq. ([7]). For a general cycle 
(B(X), (p(X)) also £l xy (k) depends on A through B (see 
the supplement). Expression (fTT|) is the general formula 
for the adiabatically pumped current when magnetic field 
undergoes a cycle in y — z plane. To be more concrete, 
we consider a cycle where the magnitude of the magnetic 
field is kept constant B(X) — B but the angle makes a full 
rotation (p(Xf)-ip(Xi) = 2ir. Then fl xy (k) is independent 
of A and the pumped charge in the cycle becomes 



' B f dk 
Q = 2ne[- / — Q xy (k) 



(12) 



The sign of Q is inverted when the field executes the 
loop in the reversed direction. Pumping formulas (|11[) . 
(fT2| remain valid in the vicinity of the topological phase 
transition despite the gap E\{k = 0) — i?-i(fc = 0) may 
become arbitrarily small because a field-induced mixing 
of the E_i(k) and Ei(k) bands is prohibited (as pointed 
out in the supplement). Validity of pumping results 
requires executing the cycle in time T which satisfies 
T- 1 < mm k (E 2 {k) - Ei(k))/H. The results in Fig.[3]a) 
show that for the chosen parameter values the pumped 
charge per cycle is of the order of Q ~ 0.5 e in the topo- 
logical phase. 

It was recently demonstrated experimentally that the 
direction of the spin-orbit field can be resolved quite ac- 
curately and indeed is perpendicular to the wire [l], 0| ■ 
Considering that in these experiments magnetic field was 
also rotated in the plane perpendicular to the spin-orbit 
field, central ingredients to test prediction (fl2|) experi- 
mentally seems to be realized. Finite size effects and dis- 
order could modify the quantitative results, but in the 
light of recent developments an experimental realization 
of adiabatic pumping (|12[) seems promising. 

Conclusion- In this paper we predicted two distinct 
magnetoelectric effects in superconducting nanowires 
with Rashba coupling, giving rise to non-dissipative cur- 
rents as a response to applied magnetic fields. The two 
effects, distinguished by whether the applied field has a 
parallel component with the spin-orbit field or not, re- 
sult from the interplay of the Rashba coupling, magnetic 
fields and superconductivity. We proposed that the pre- 
dicted effects could be observed in recently realized ex- 
perimental setups and employed in probing the topologi- 
cal phase transition and realizing adiabatic charge pump- 
ing. 

It is pleasure to thank Jay Sau and Takuya Kitagawa 
for discussions. The author acknowledges Academy of 
Finland for support. 
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SUPPLEMENTARY MATERIAL 



Here we present supplementary results and formulas that were cited in the main text. 



Magnetoelectric effect for parallel fields B x 7^ 



In the main text it was noted that application of magnetic fields with a non-vanishing B x component (which is, 
according to our conventions, parallel to the effective magnetic field arising from the spin-orbit coupling) give rise to 
current 



e f dk \ / hk a , , 
J =2 ^ - + ft Tr ^>/^o + ^ 



(13) 



where Pi = Ylj^i e -e- are ^ ne projection operators to the four eigenstates Ei of the Hamiltonian (JTJ) of the main 
text and the trace should be calculated over the Nambu and spin indices. Current components Jo and J\ correspond 
to the first and the second term in the curly brackets in Eq. (fT3")) . Due to the finite parallel component B x ^ 0, the 
energy bands do not admit analytical expressions. The explicit formulas for the traces in terms of Ei are given by 



Tr [Pi<j x ] = 



so Ji becomes 



Ji 



ea 
2h 



dk 



Yl&iiEi-Ei) 



4/, 



2;„2 



Bt + 2e k akEi + B x 3A Z + 3ei + 3a z k 



B 




^^rx[^ i {Ei-E j 



Bt+B x \ 3A^ 



3a 2 k 2 



B + Ef 



2 ^ 

k 



Ei 



2ehakEi 



(14) 



(15) 



From the property Ei(k,B x ) — Ei(—k,—B x ), it follows that current (fl~3|) is an odd function of the parallel field 
J(B X ) = —J(—B x ) (implying that J(0) = 0). Calculation of current generally requires solving the energy bands 
numerically and then evaluating Jo and J\. However, to evaluate current in the linear order in B x one can make 
further progress analytically. By applying the first order perturbation theory, the energy bands can be solved as 
E±i(B x ) = ±E X {B X =0)-SE and E± 2 (B X ) = ±E 2 {B X = 0) + SE where 8E = ^.fffi and E % (B X = 0) are given 
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by Eq. @ in the main text. Current up to the linear order in B x can then be evaluated by inserting the expressions 
for energies in (j 1 5|) and retaining only linear terms. This procedure yields 

eaB x fdk^ 4ft (-1)' 
_ eaB x f dk 4(l-2/Q(-l 

The spin-independent part Jo does not contribute to current in the linear order in B x so J = J\. Defining D = 
E\ — E\ = Ay / B 2 (el + A 2 ) + e'f.a 2 k 2 , Eq. (fT())) coincides with the current formula (g} in the main text. 



(-iy 



4E? 
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:E 2 
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El 
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(16) 



Magnetoelectric effect for perpendicular fields and adiabatic pumping 

As discussed in the main text, the perpendicular magnetoelectric effect (for fields with B x — 0) and adiabatic 
pumping is determined by the quantity 

■w«-(-t)i:t^*-w (it) 

x 7 n,n' 

where P n is a projection operator to the eigenstate of the Bogoliubov- de Gennes Hamiltonian with eigenvalue E n 
given by Eq. ([2]) in the main text and /„ is the Fermi function at energy E n . The projection operators are given by 

Ei = rij^i E -E- wn i cn m this case simplify to P±i = ^ ( 1 ± -^-J E i_ E l and analogously for P±2 by interchanging 
subscripts 1 O 2. Evaluation of the traces in Eq. (jXTJ) is now in principle straightforward but a tedious task. The 
explicit form of £l xy is given by 

( A 2 + £ 2 fc E l+ E 2 \ (h + f 2 -l) ( A 2 +e 2 E 2 -EA (A - f 2 ) 1 , 
\ E x +E 2 4 J {E x +E 2 ) 2 \E 2 -E l 4 J {E 2 -E l ) 2 \ ' 1 J 

This form shows that VL xy vanishes if either B or a vanishes. In the linear response current (Eq. ([8]) in the main text) 
B plays the role of the static part of the field while in the adiabatic pumping process (Eq. f| 12[) in the main text) 
B = B(X) is the slowly varying function of the adiabatic parameter. In the adiabatic pumping process the energies 
depend on the adiabatic parameter through B(X). 

It is worth noting that in Eq. (|17p the terms in the sum corresponding to pairs n — 1/2, n' = —1/ — 2 vanish, 
since the system possess a selection rule (E n (k)\<r y / z \E_ n (k)) — 0. This selection rule can be shown with the help of 
the chiral symmetry {H,a x T y } — (which is present when B x = 0). Therefore perpendicular magnetic fields cannot 
couple the E\ and bands and VL xy remains analytic at the phase transition where E-\{k) = E\(k) — at k — 0. 
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